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Abstract 

Generation of sound due to scrubbing of a jet flow past a nearby solid surface is 
investigated within the framework of the generalized acoustic analogy theory. The analysis 
applies to the boundary layer noise generated at and near a wall, and excludes the scattered 
noise component that is produced at the leading or the trailing edge. While compressibility 
effects are relatively unimportant at very low Mach numbers, frictional heat generation and 
thermal gradient normal to the surface could play important roles in generation and 
propagation of sound in high speed jets of practical interest. A general expression is given 
for the spectral density of the far-field sound as governed by the variable density Pridmore- 
Brown equation. The propagation Green’s function should be solved numerically starting 
with the boundary conditions on the surface and subject to specified mean velocity and 
temperature profiles between the surface and the observer. The equivalent sources of 
aerodynamic sound are associated with non-linear momentum flux and enthalpy flux terms 
that appear in the linearized Navier-Stokes equations. These multi-pole sources should be 
modeled and evaluated with input from a Reynolds-Averaged Navier-Stokes (RANS) 
solver with an appropriate turbulence model. 


1.0 Introduction 

The Fundamental Aeronautics program at NASA faces new technological challenges in order to 
meet recent air transportation rules that include reduced noise, emission, and fuel consumption. 
The U.S. commercial aviation is projected to increase by 70% in the period 2010 - 2030, 
resulting in more exposure to air traffic noise in the communities surrounding airports. Recent 
noise metrics require a cumulative (cutback, sideline, and approach) noise reduction in the 
subsonic transport of 32dB by 2015 and 52dB by 2025 (referenced to 737-800 aircraft with 
CFM56-7B engines). In 2008, NASA instigated a call for conceptual design of future 
commercial transport in order to meet specific goals related to noise, air pollution, and fuel 
consumption. Early studies indicated that both propulsion system and engine placement have to 
be modified in order to achieve these new objectives. In particular, placement of the aircraft 
engine for the purpose of reduced community noise poses new challenges in structural design, 
material selection as well as computations. Configurations such as over the wing engine mount 
(Fig. 1), distributed propulsion, Hybrid Wing Body (HWB) concept, and/or high aspect-ratio 
rectangular exhaust geometry with extended beveled surfaces designed to shield noise from 
reaching the ground, also produce new noise sources due to scrubbing of the flow past structure, 
skin vibration, boundary layer noise, and edge noise due to the scattering of sound from sharp 
edges. 

In 2010, a series of tests were initiated at NASA Glenn Research Center in order to study the 
propulsion/airframe integration. These experiements, under the umbrella of “Jet-Surface 
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Interaction Test (JSIT),” were intended to guide analytical studies towards prediction model 
development. The initial tests consisted of placement of a flat plate in the proximity of a round 
jet. Details of the experimental setup and test cofigurations are provided by C. Brown [Ref. 1], 
and G. Podboy [Ref. 2], Measurements were gathered on both sides of the plate, commonly 
referred to as the reflected and shielded sides. In addition to the far-field microphone noise 
measurements [1], phased array source localization experiments [2] were carried out in order to 
shed light on the noise source location and its frequency content under different exhaust 
conditions. 

An example of JSIT noise data that shows 
power specral density per Strouhal frequency 
in a 2-in diameter (D =5. 08cm) Mach 0.90 
unheated jet, as reported by Brown in Ref. 1, 
is reproduced in figure 2. Here a solid surface 
extending 12-diameters from jet exit in the 
downstream direction was positioned at DU 
from the jet centerline. Far-field measurements 
are shown at 90° and at 1 00D for a jet in 
isolation (no solid surface), as well as shielded 
and reflected sides of the jet in the presence of 
the surface. 

An extensive database has been generated as a function of the jet operating conditions and the 
length and position of the nearby solid surface. Similarly, an example of source localization data 
for the above jet is reported by Podboy as figure 17 of Ref. 2 (not shown here). It presents a 
beam-form map of the source location for an isolated jet as well as when surfaces of varying 

It is commonly accepted that the 
trailing edge noise is responsible for 
the enhanced low-frequency and 
peak noise amplitude measured on 
both reflected and shielded sides of 
the plate as seen in figure 2. A 
distinctive feature of the trailing 
edge noise component is its cardioid 
directivity shape and a dependence 
on jet velocity to the fifth power 
[Refs. 3, 4], Recent measurements 
by Lawrence et al. [Ref. 5] indicate 
that the velocity power factor could 
be as large as 6 at longer plate 
length L that yield a wetted trailing 
edge. FFowcs Williams and Hall 
[Ref. 3] describe a solution to the 
edge noise from a rigid plane in the 
presence of a turbulent flow using 
Lighthill’s acoustic analogy. They 


length are placed between jet and the phased array. 



Figure 2. Sample JSIT noise measurement as 
reported by C. Brown in Ref. 1. 


4 * 



Figure 1. Jet surface interaction. 
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consider a very low Mach number flow, and take advantage of an appropriate Green’s function 
applicable to a point source in the vicinity of a half-plane [Ref. 6]. Howe [Ref. 7] reviews 
various theoretical approaches in understanding the trailing-edge noise, and summarizes diverse 
theories into a single model that exhibits the dependence of the acoustic field on the turbulent 
fluctuations near the edge as well as the significance of the Kutta condition at this point. 
Goldstein [Ref. 8] includes the effect of the mean velocity gradient on the trailing edge noise by 
imposing a convecting vortex source, or gust, of the general form G(x, /U(x 3 ) — t,x 2 ,x 3 ) on an 
incompressible transversely sheared mean flow. 

Our interest here is the turbulence-generated noise due to the scrubbing of the jet flow past an 
adjacent flat surface. At very low Mach number, Howe [Ref. 9] describes the far field acoustic 
domain in terms of the wall-pressure wave number- frequency spectrum. Various models have 
been suggested for wall-pressure spectrum [Ref. 10] that are used in conjunction with a 2D or 3D 
free-space Green’s function to express the far-field spectral density. These approximations 
ignore the mean flow refraction, compressibility effects, and sound produced outside the 
Boundary Layer (BL) at distances that are not small relative to the acoustic wavelength. At jet 
velocities of practical interest typical of a jet exhaust, frictional heat generation and thermal 
gradients in the BL could play important roles in generation as well as propagation of sound. 
Additionally, skin vibration could influence the surface boundary conditions and the subsequent 
radiated sound. 

In section 2 we derive an expression for the propagation Green’s function (GF) applicable to the 
scrubbing noise on a flat surface. The GF needs to be obtained numerically when the mean flow 
profiles are non-uniform. An expression for the far-field pressure is given in section 4. Ideally a 
RANS solution for the flow in the proximity of the surface provides the mean flow and 
turbulence information required in evaluating both the source strength and the GF. Sample GF 
computations using analytical profiles for the mean flow and a summary are provided in sections 
6 and 7. 


2.0 Formulation of the Scrubbing Noise 

Consider turbulence-generated noise due to the scrubbing of a jet flow with an adjacent flat 
plane (Fig. 3). The mean flow is represented as a unidirectional transversely sheared flow 
(independent of xi and x 2 directions). We consider the mean static pressure p within the BL to 
be a constant. At very low Mach numbers the compressibility effects may be negligible, 
however at Mach numbers typical of a jet exhaust, frictional heat generation and thermal 
gradients in the BL could play important roles in generation as well as propagation of sound. 
The RANS solution provides such effects via the action of viscosity near the surface. Hence the 
noise generation process due to the scrubbing of a jet with a nearby surface is, in general, 
described by full acoustic analogy as considered in Appendix-A, and by Pridmore-Brown 
equation within a locally parallel flow approximation [Ref. 1 1] 


Ln' = T, n' = 
where L is the operator 


p'(x,t) 

yp 


(i) 
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L = D 


d 2 -JT (c2 JT 

dXj dXj , 


| 2 du a 2 D= d_ +U _<L 

dx . 5xj dx . ’ 5? dx t 


( 2 ) 


and the source term r(x,f) on the right hand side of (1) is defined according to the generalized 
acoustic analogy [Ref. 12] 


r = 


yp 


d , de,., , 

-D—(c 2 -^) + 2c 2 
ox , ox. 


G = -(y-i) 


2 


dx. 


dXj dx, dx. J 

(3a) 

riTI ^ 

(/ pv'h' 0 ) + {pv[v ')— , 

(3b) 


stress tensor e tj is the difference between the fluctuating and Favre-averaged Reynolds stresses 


e ij = -ipv'v'j - pv'v'j) , (3c) 

and h ' 0 denotes the fluctuations in the moving- frame stagnation enthalpy (see the nomenclature). 



Figure 3. Jet scrubbing on a nearby surface. 


The mean velocity U and density p are functions of the normal coordinate (X 3 >0), and c is the 
mass-averaged speed of sound 

U = U(x 3 ), c 2 =c 2 (x 3 ) = yp/p(x 3 ). (4) 

Note that, in practice, both the mean flow and turbulence depend on the transverse direction xi as 
the source strength eventually diminishes at some distance ±L 2 away from the jet centerline. In 
the stream-wise direction, the scrubbing effect may be limited to a distance 0 < x 1 < L { , or the jet 
may persist beyond the plate edge at L\ where interaction of the turbulence with the trailing edge 
produces edge noise. Here, in formulating the scrubbing noise problem, we limit the source to a 
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finite X 2 span by considering the mean flow as locally parallel in both xi and X 2 directions - thus 
express the problem in three dimensions. This assumption, although more suited to jets with a 
high aspect-ratio rectangular exhaust, permits source-volume integration within the wetted span 
in X 2 direction. A round jet may be transformed into a rectangular strip by using an appropriate 
conformal mapping transformation. When all changes in span-wise X 2 direction are neglected, the 
sound intensity should be expressed as per unit thickness. The special 2D case will be discussed 
in section 5. 


We seek a solution to equation (1) of the form 

K'(x,t) = j jG(x,t;y,z)r(y,z)dzdy , (5) 

y f 

where G denotes the GF 

LG(x,f,y,z) = S(x-y)8(t- r) . ( 6 ) 

In a bounded media, the above form of the solution requires for the GF to satisfy certain 
boundary conditions on the plate. The source-volume integration in (5) is carried out over noise- 
generating region y in the vicinity of the scrubbing surface and at the emission time T . The 
solution described in the following excludes the trailing edge noise that travels to the opposite 
side of the surface as alluded to earlier. This latter noise component requires a Wiener-Hopf 
type solution. 

We define a Fourier Transform (FT) pair with respect to (x p x 2 ) coordinates as well as time t. 
The corresponding wave number vector is denoted as k t = (k t ,k 2 ) 


G(k t ,x 3 ,y 3 ,co) = ,t\y,T)e ' 


[ k x {x x -y ] )+k 2 (x 2 -y 2 )-co(t- T )] 


d (X| — V| )d(x 2 -y 2 )d(t-T) 


G(x,t,y, T)= J J J G(k t ,x v y 3 ,w)e ilk ' ( - Xl ~ yi)+k2l - X2 ~ y2) ~ 0>(t ~ r)] dk l dk 2 d( 0 . 


( 7 ) 


Throughout, an over hat carrot is used to denote a FT variable as defined above. 


The transform of equation ( 6 ) is 


dxl 


( /ji 


- + 


(c 2 )' 2k x U' 


c -co + k x U 


dG 

dx. 


- + 


( CO + kJJ) ,2 j2 


2 ^ 

— k , 2 — i 

J 


G • ' ! . ( 8 ) 


(Inf c 2 (—(Q + kJJ) 


where prime as used on mean variables c 2 and U denotes differentiation with respect to X 3 . The 
equation possesses a critical layer atx 3 where factor -co + k t U (x 3 ) vanishes. To avoid this 
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singularity we require k { <(0 /U e where U e denotes some peak mean velocity in the BL. This 
upper limit for k\ is usually referred to as the non-convective domain for the wave number. As 
discussed in [Refs. 13, 14], in the convective domain k x >0)!U e it becomes necessary to provide 
a loss mechanism, such as turbulent viscosity, in order to evade the singularity. The dominant 
value of k A contributing to the far-field noise will be determined in section 3. It will be shown 
that when U e < the upper limit condition stated for k { will be satisfied at all observer location. 


For convenience, the second-order linear differential equation (8) is rearranged as 


+ f(k t ,x 3 ,co)G = 8(x 3 -y 3 ) . (9) 

Since equation (9) is self-adjoint, the newly defined GF, G , is symmetric in (X3, yj), i.e., 
G(k t ,x 3 ,y 3 ,0 ) ) = G(k t ,y 3 ,x 3 ,0) ) and is related to the GF of interest G as 


G(k t ,x 3 ,y 3 ,co ) 


i 1 -co + kJJix^ 
(2 nf c(y 3 )c(x 3 ) (-m + k 1 t/(y 3 )) 2 


G(k t ,x 3 ,y 3 ,co ) , 


( 10 ) 


where 


f(k t ,x 3 ,( 0 ) = x 2 ~ C 2 - 


(-co + kJJ) 2 


+ 


kjr_ | (c 2 y kjr_ 

-to + kjC/ c 2 -to + kjf/ 


X 2 (k t ,x 3 ,co) = - — k\-k\, 

c 


1 d 

f(c 2 Y ) 

l 

+ — 

f (c2) 'l 

2 dx 3 

v c2 J 

4 

v c2 J 


( 11 ) 


and (c 2 )'/c 2 = (dp /dx 3 )/ p . 

It should be noted that unlike the present problem, the propagation equation might not, in 
general, reduce to a self-adjoint from. Such is the case when the governing equations are written 
for a spreading flow [Ref. 15]. In that case, working with the adjoint GF will entail additional 
integrals that should be carried out on the surfaces. A general solution for the five field variables 
using a vector GF formulation is provided in Appendix A. The second integral in equation 
(A16) in the appendix describes the contribution to the acoustic field due to the nearby surfaces. 
In the formulation described herein, the GF satisfies conditions that make the surface integral 
equal to zero. 

Far away from the noise generating region, as x 3 — » °° the mean velocity and the sound speed 
approach their respective ambient values U m and cy , and equations (9) and (11) show that 


3 2 G 

d 2 x; 


+xiG= 0, 


x 3 ->°° 


(12a) 
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where 


xl = (-K 0 + k\Mj 2 -k;-k;- k 0 = cq/c„, M„ = ujc„ (12b) 

The Mach number takes account of the flight effect when surface is in motion relative to the 
ambient. The far-field solution to equation (12a) is 

G(k t ,x 3 ,y 3 ,co) = b e~' x ~ x \ x 3 — (13a) 


where complex number b has dimension of length, and in general depends on (0,k t ,k 2 . An 
outgoing propagating wave requires the negative root of when xt > 0 , and when xt is 
negative the branch-cut is chosen such that x~ = - *0 xt 0 12 • When xt < 0 the pressure decays 
exponentially normal to the surface (evanescent waves). Aerodynamic noise scattered off the 
edge of a half-plane is sometimes described [Ref. 16, 17] as the interaction of evanescent waves 
(rather than propagating acoustic waves) with the edge. These disturbances are dominated by 
spectral components that possess a subsonic phase velocity in a constant X3 plane. 


Since we are interested in the scrubbing noise that reaches a distant observer, the radiation field 
would be dominated by wave numbers, in constant X3 plane, that possess a supersonic phase 
velocity 


co 


M 




+ k: 


> C 


(13b) 


The far-field boundary condition for a bounded solution to equation (9) requires 


dG 

dx. 


+ iXJG = 0, x 3 ^ °° , 


(14) 


and a no slip boundary condition on the scrubbing surface (U = 0, x 3 — » 0) combined with the X3 
component of the linearized momentum equation within the turbulent boundary layer (TBL) 
shows that 




+ - 


= 0 , 


^ 0 


(15) 


dt dx 3 

Upon replacing p with the mean fluid density p o near the surface, we apply a FT to the latter 
equation to obtain 


O dp' 

-i«)p o v 3 + — ^0, 

ox. 


x 3 — ^ 0 


(16) 


Assuming that Fourier components of the normal velocity and pressure on the surface relate 
through the surface impedance Z(co) as 
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(17) 


p' = Zv' 3 , x 3 = 0 


we find 

dn' 

dx 3 



*3=0 • 


(18) 


Since the FT of a convolution is the product of the Fourier transforms, equations (5) implies 


n'(k t ,x 3 ,co) = (27t) 3 J G(k t ,x 3 ,y 3 , (o)T{k t ,y 3 ,(o)dy 3 , 

ft 


(19) 


and upon using (19) in (18) we find 


A.!M.g = 0 

dx 3 Z 


* 3 =° 


( 20 ) 


Our interest is in the boundary condition applicable to G rather than G . Substituting equation 
( 10 ) in ( 20 ) shows that 


y/G - 0 , -x 3 =0 

OX 3 


yf(k l ,a>,Z) = 


!S- 

Z c; 


c' k 

+ ^ + ^U'( 0 ) , 
c 0 O) ) 


Z(co) = 


Z((0) 

P C 

§ oo oo 


( 21 a) 

(21b) 


where c 0 is the mean sound speed in the fluid at the surface interface, and c' a is its derivative at 
X 3 = 0. An alternative way of deriving the surface boundary condition (21a) is described in 
Appendix B. 

The problem at hand is thus reduced to solving equation (9) subject to the two boundary 
conditions (14) and ( 21 ). 


The two boundary conditions could be expressed in a more general form as 


3 k ' -> 

+A* =Yi > 

ox. 


i - 1,2 


( 22 ) 


where i=\ and 2 refer to the boundary points at X 3 = 0 and 00 . A non-zero value of y t (such as 
presence of a source term in the linearized momentum equation near the surface) would 
contribute to the homogeneous solution to the differential equation. However, for our problem 
we consider y l = y 2 = 0 , thus the homogeneous solution is identically zero, and we are left with 
the particular solution only. 
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The GF in equation (9) may be expressed in terms of two linearly independent solutions 
V i (k t ,x 3 ,co ) to the homogeneous equation 


d 2 V i /dx 2 3 +fV i = 0, i = l,2 


(23a) 


such that V\ and Vi satisfy the homogeneous boundary conditions at 0 and oo, respectively. For 
brevity we omit k t and co from the arguments of V\ and Vi. 


and 


Mg) 

dx 3 


-i//V 1 (x 3 )=0, x 3 = 0 


M* 3 ) 

dx 3 




X 3 — > oo 


(23b) 

(23c) 


G(k t ,x 3 ,y 3 ,co) = 


V 2 (jc 3 ) V, (y 3 ) / W (y 3 ) , y 3 <x 3 
yi(x 3 )V 2 (y 3 )/W(y 3 ), y 3 >x 3 


(24) 


and IF (V, ,V 2 (y 3 )) is the Wronskian. For convenience we normalize both V\ and Vi at zero 
argument y 3 = 0 such that 


W) = V 2 (0) = 1, 


(25) 


and since according to Abel’s formula the Wronskian to equation (23a) is independent of y 3 , 
then we evaluate W = VF at the surface. Equations (23b, c) and (25) conclude that 


W o (k t ,CQ,Z) = v;(0)-xi/(k l ,0),Z). (26) 

The above expression shows that the Wronskian W 0 depends on the BL velocity and temperature 
profiles aty 3 = 0, as well as the surface impedance Z through parameter y / . It is noted that 
equation (23 a) supports Helmholtz instabilities that would arise when k t and co satisfy the 
Eigen-wave condition W a (k t ,co ,Z) = 0 . The instability waves could be triggered when the BL 
profile supports the Eigen-wave condition. These waves could grow exponentially and dominate 
the region they occupy. Subsequently we need to set the wave number limit to avoid these 
special cases. For example, when the surface is rigid and the BL is considered as incompressible 
such that ii/(k l ,co,Z) = k l U'(0)/ co , we should require V '_(()) I U'(0) >1/U e for the Wronskian W 0 
to remain non-zero for all wave numbers within the non-convective range k x <(0 IU e . However 
as we will discuss shortly, F 2 '(0) is not known a priori and is determined only after equation (20a) 
is solved as a boundary value problem. Jones [Ref. 18] studied the instability of a two- 
dimensional shear layer using a model velocity profile in which the velocity increased linearly 
from 0 to U e over a distance of h from the surface and then remained constant. He showed that 
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the onset of Helmholtz instabilities occur at a Strouhal frequency of coh / U e ~ 1 / 3 where co 
denotes the frequency of sound waves falling on the shear layer. At frequencies above this 
critical value the layer remained stable. 


Upon placing (24) into (10), the Fourier transformed GF at an observer point above the source 
region is given as 


G(k t ,x 3 ,y 3 ,a > ) 


i 1 -O) + k A U(x 3 ) V 2 (x 3 )V t (y 3 ) 

(2 nf c{y 3 )c{x 3 ) {-(o + k l U{y 3 )) 1 W a 


x 3 > y 3 


(27) 


An alternative form of the GF valid for x 3 <y 3 (see equation 24) should be used if one were 
interested in the surface pressure. The two functions V i are determined numerically. The first 
variable is evaluated when equation (23a) is solved as an initial value problem with 1^(0) = 1 and 
1^(0) = y / . The second variable is determined when (23a) is solved as a boundary value problem 
subject to F 2 (0) = 1 , and V 2 (x 3 ) + ( v 3 ) = 0 as x 3 — > °° . Fortunately, as we shall see next, 

there are only certain values of wave numbers k t and k 2 that make the major contribution to the 
radiated far-field noise. 


3.0 Stationary Phase Method 

We apply an inverse FT with respect to wave numbers k { and k 2 and write the GF as 


G(3c,y;©)= J J G(k, ,x 3 ,y 3 ,co ) e' lk ' (x '- y ' ,+k > ix ^ ,l dkylk 2 


(28) 


Since the interest is in the far-field noise, the appropriate from for V 2 is 

V 2 (k t ,x 3 ,0J) = b 2 e ~' x ~ X3 , x } “ . 


(29) 


Complex number b 2 (k t ,co ) has a dimension of unity, and should be a known parameter once the 
boundary value problem (23a) is solved for Vi. Placing (27) and (29) into (28), the double 
integral involving wave numbers k\ and ki becomes 


G(x,y;«) = 


-(U + k 1 t/(x 3 ) b 2 V 1 (k t ,y 3 ,co ) 
(2nf d'v irf.Y.) J J (-oj + t/Mvjr H'.it.oi.Z) ‘ 




t iO(k t ,x ,co) 


dk ] dk 2 , (30) 


where 


e(k t ,x,co) = k l (x l -y 1 ) + k 2 (x 2 -y 2 )-x^x 3 . (31) 

We adopt a spherical coordinate system where 6 is the angle the observer position vector makes 
with the direction of the mean flow (Fig. 4) 
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(jtj -yi,x 2 -y 2 ,x 3 ) = R(sin</>cos0, cost/), sin(/)sin6), x 3 — , 


(32) 



Figure 4. Spherical coordinates 0 <6 <7t, 0 <(j)<jr . 


where R = \x - y \ . The wave number vector in expressed in a similar fashion as 

= K 0 (sm(j)' cosO' , cos0', sin^'sin#') . (33) 

Using (32) and (33) in (31), and upon neglecting the flight effect, the phase factor becomes 

0(^,x,ft)) = i?X' o (sin^sin0'cos(0-0') + cos^cos^'), M m = 0, (34) 

where we have selected =-K 0 sm(j) , sm6' . The large parameter in applying the stationary 
phase technique [Ref. 19] is Rk o »1 . This requires the source to observer distance R to be 
much larger than the wavelength of the radiated sound. If we define a Strouhal frequency as 
St - co8 0 /U e , where 8 a is the BL thickness (distance from the wall at which the mean velocity is 

0.99 U e ), then the large parameter requirement implies St » c oo S 0 / (U e R) which is not a very 
stringent condition even at relatively low frequency. 

The Jacobian in making coordinates transformation from (k p k 2 ) to ( Q\(j )' ) is K^sin0'sin 2 0' or 
dk { dk 2 = k 2 0 s\n 6' sm 2 <j)'d0'd<j)' . 

Using the new coordinates, the point of stationary phase (6* ,(f) ) is obtained when we set equal 
to zero the two partial derivatives of phase factor 0 . 


30 ^ 30 ^ 

(35) 

— = 0, — = 0 . 
33' 30' 

which results in 


II 

V- 

Cb 

II 

Cb 

(36) 
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Equation (36) simply implies that in the absence of flight effect, the major contribution to the 
radiated sound field is due to wave numbers with magnitude k 0 in the same direction as the 
observer point (#,</>) . Subsequently the contribution to the integral in (30) in the vicinity of the 
stationary point is obtained when we expand the phase 0 (to second order) about this point, and 
evaluate the rest of the integrand at the point of stationary phase k' = (k’,k 2 ) . 


k* = x o (sin 0 cos#, cosrp). 


G(x,y;m) 


i -to b 2 (k' , co ) V, ( k' , y 3 , co ) ^ 

(2^) 3 (-to + yt;t/(y 3 )) 2 c^c{y 3 )W 0 (k[,co,z) 


JJ ex P 

a' Ai' 



l-^(0'-0) 2 sin>-|(f-0) 2 


(k 2 0 sin Osin 2 ^)dO'd(j)'), 


(37) 


as (R — > <=o) . The above double integral is now evaluated as 


e lK ° R K 2 0 sin Osin 2 0 J J exp - ( a 2 sin 2 <j) + /5 2 ) dadj5 - -2niK a sint/sin ^ e K « R ( 38 ) 


sin 0 sin t 


R 


Placing this result back in (37) shows 

1 e ,K ° R sin 0 sirup b 2 (k* ,co)\ j (k' t ,y 3 ,(0 ) 


G(x,y;to) — 


1 


(2tt) 2 /? clc(y 3 ) W a (k;,co,z) 


f tt / x \ 2 ' 

l-sin0cosO 


,^(y 3 ) 


(39) 


V 


J 


The GF G(x,t;y , T ) is now obtained as 


G(x,t;y, t)= J G(x,y;to)e ,a,( ' T) dto . (40) 

Note that with A", = k o sin^cosO , the upper limit placed on k A in section 2 will be satisfied at all 
subsonic conditions U j I c rx < 1 . At supersonic conditions only certain observer locations will be 
subject to the critical layer singularity. 


4.0 Far-Field Acoustics 

The far- field sound spectral density is evaluated as a FT of the pressure auto-covariance 


°° r i 

p 2 (x,co) = (yp) 2 J e ,on — J K\x,t)n\x,t + T)dt 


\dr 


(41) 


where T denotes some large time. Upon placing (5) and (40) into (41), it is shown that [Ref. 20] 
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( 42 ) 


p 2 (x,co)=jj J G\x,y-ll 2 ,( 0 )G{x,y+ll 2 ,( 0 )q{y,l,T)e iw dTdldy , 

y | T=^= 

where q denotes a two-point space-time correlation between noise generating sources at points A 
and B separated by spaced and timer . Using the source term in equation (1) we have 

1 T 

q(y,lz)^(Yp) 2 lim—\T*(y-l/2q)ny+l/2,t + T)dt . (43) 

T ^°°2 T J T 

Since the exponential phase factor in the GF equation (39) depends on the distance R=lx-yl 
between observer point x and source point y , the product of the GF and its conjugate (denoted 
by *) is evaluated at the center of the correlation multiplied by a phase factor exp (-ik.%), where 
wave number k is directed as x—y and has a magnitude of k 0 


p 2 (x,co) - 1 1 G(x,y;«) I 2 J j q(y,^,T)e im - ik l dzd^dy . 

y l t=-°° 


(44) 


As is a common practice in acoustic analogy noise-prediction models, the space-time FT of the 
source correlation function, denoted as the inner double integral in (44), is evaluated in closed- 
form once an acceptable physics-based source model is obtained. The modeling becomes much 
less cumbersome if the spatial derivatives present in the source were to be transferred to the GF 
within equation (5) and prior to forming the sound spectral density function. Subsequently, the 
final GF as defined in (39) could be subject to new spatial derivatives at the source point y while 
source T(y,t) would be free of similar derivatives. 


Since the GF depends on both V\ and Vi, which in turn depend on co ,0 , and (f ) , the pair ( V t ,V 2 ) 
is determined numerically per pre-selected observer angles and frequency using the local BL 
mean velocity and density profiles. The integral over the source region in equation (44) is then 
evaluated to obtain the far-field noise as a contribution from independent correlation volume 
elements that comprise the source. As pointed out earlier, the volume integration is limited to 
regions with relatively stronger source strength while the mean flow is considered as locally 
parallel in both xi and X 2 directions. 

In the limit of very low Mach number (M <0.10), it is argued [Ref. 9] that there is a direct 
correspondence between the acoustic domain of the wall-pressure spectrum and the far-field 
spectral density of the radiated sound. This approximation ignores mean flow refraction and 
sound produced by the turbulence outside the boundary layer at distances that are not small 
relative to the acoustic wavelength. Following this argument, at large distance R from a hard 
wall region of area A, Howe [Ref. 21] proposed an expression for the sound spectral density in 
terms of the so-called blocked pressure, which in the present notation is 


P 2 (x,cq) = Ak 2 0 


sin 2 6 sin 2 (j) 

R 2 


P(k t ,a>) , 


(45) 
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where P(k t ,co) is the wall-pressure wave number-frequency spectrum, and is defined as a FT of 
the space-time pressure correlation function < 7 (£j,£ 2 ,T)of the wall-pressure 

p (k t ,Q)) = -—^ J J q(^ ,T)exp(/mr - ik t .l)dzd^d^_ . (46) 

Here is the spatial separation vector of the correlation on the surface and k t is the 
corresponding wave-number vector. The wall-pressure correlation function over a rectangular 
area A is defined as 


^ ^ T L x L 2 

^i,4,T)=lim — - J J j p'(y l ,y 2 ,0,t)p'(y l + ^ l ,y 2 + ^ 2 ,0,t + z)dy l dy 2 dt . (47) 

Various models for P(k t ,(o) have been proposed in the literature, however there is no universally 
accepted form for the wall-pressure spectrum that could directly be substituted into (45). For 
example low wave number measurements of Sevik [Ref. 10] suggest 


P{k t ,(Q) niWJcJiyJU') 

plvlSl (c o8jU e r 


24 <^<240, 0.01< — <0.15. 
U. c 


(48) 


where v* ~ 0.03£/, is the friction velocity, and S a denotes the BL thickness as defined earlier. 


More complicated expressions have been suggested for the wall-pressure wave number- 
frequency spectrum P(k t , co) with additional parameters that are meant to account for the 
attenuation and refraction of sound in the BL [Ref. 22] and/or surface roughness [Ref. 23] and 
dependence on the individual wave number components on the surface. The validity of these 
empirical models remains to be confirmed experimentally. A comparison of semi-empirical 
equation (45) with the analytical solution given in equation (44) shows that the latter should 
account for all of the above-mentioned effects when source information as well as the GF are 
available within and outside the BL. Ideally, a RANS solution to the wall-bounded flow would 
provide both the mean flow and turbulence information required in evaluating the GF and the 
source correlation in a turbulent boundary layer and the main jet. 


5.0 Two-Dimensional Approximation 

When mean flow and turbulence are independent of the transverse direction (unlike in 
section 3 where the rectangular flow was limited to a finite X 2 domain), the analysis is carried out 
in two dimensions. Following the steps outlined earlier, the new GF is obtained from an 
equation similar to (30) 


G(x,y;m) 


[ \ r -gj + kpjxs) bjVj (k t ,y 3 ,m) j ea , 

(2n) 2 c(y 3 )c(x 3 )“ ^—co + k l U(y 3 )) 2 W 0 {k„( 0 ,Z) 


(49) 
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where phase parameter is now defined as 


Q(k l ,x,co)=k l (x l -y l )- x^x 3 , xl = K l ~ k i • 

It is readily shown that the point of stationary phase is k[ = k 0 cos 6 , and 


( 50 ) 


G(x,y;«) — 


(2n) m ^JK^R clc(y 3 ) W 0 (k[,co,z) (1 _ ' 


(51) 


The source correlation now uses | = (|,,£ 3 ) as a separation vector, and the far- field sound is 
evaluated using equation (42) with the volume element as dy = dy l dy 3 , and per unit flow 
thickness in X 2 direction. Note that both 3D and 2D expressions for the GF could be written in a 
similar from if equations (39) and (51) were normalized with respect to their corresponding free- 
space GF. 


6.0 Numerical Results 

Sample Green’s function calculations are presented using analytical representations of the jet 
mean velocity and temperature profiles in the proximity of a flat surface. 

The mean velocity profile adopted here depends on the normal direction y 3 only, and is modeled 
for a jet with thickness ZT - 2 in 


U(rp 

u j 


tanh( 


M), 


11 u * 

— + — tanh — 
2 2 


( 


V 


1/2 77-1 

77-I 1/2 


77 <1.05 
77 >1.05 


(52) 


Using parameters (d t ,d 2 ) = (0.10, 2.0) , the velocity profile, divided by jet velocity U ] , is plotted 
(Fig. 5) as a function of the non-dimensional distance from the wall 77 = y 3 l D r The boundary 
layer thickness for this flow is 8 a / D } - 1 .324J, . 
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Uirf) 


Uj 



Figure 5. Mean velocity profile. 


The mean static temperature is modeled as a composite of the two profiles - the first is obtained 
by placing the above velocity profile in Crocco-Busemann law 


jm 

71 


= 1 + (T r -1) 


U(V ) 


7-1 


U (77) 


V 


(53a) 


where T R is the jet plenum stagnation temperature ratio. The second profile is intended to 
simulate frictional heat generation in the proximity of the wall 


T_m 

T o o 


1 f 1 1 u 1 , 1 

— — + — tanh — ( 

d 3 y2 2 d 4 DjTI 


\ 

-DjV) . 

J 


(53b) 


Using parameters (d 3 ,d 4 ) = ( 4.0, 3.0) the temperature profile T(^) = T l (r\) + T 2 (r\) is shown in 
Fig. 6 subject to jet exit values of Uj / c^= 0.90 , and T R = 3.0 . The ideal gas law c 2 = y%T was 
used to construct the speed of sound as shown in Fig. 6. 
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Figure 6. Mean profiles for static temperature (solid line) 
and sound speed (dashed line) at T R = 3.0 . 


The GF is solved numerically at specific values of the Strouhal frequency St Q defined as 


1 cqDj 


In Uj ’ 


(54) 


We define non-dimensional parameters K -K 0 Dj-2nSt 0 (Uj / c„) , c = c / , U = U / c m , 

R- R/ Dj, and divide the 3D and 2D expressions in (39) and (5 1) by their respective free-space 

values of(-e ,KR / 4nR) and (- e lKR+,Km / VS/nc/? ). Additionally constant 1 1 (net) is factored out 
in order to define a normalized GF applicable to the far field 


G N (Ji,co) = 


sinO sincp b 2 (k* , coW^k* ,r] ,co) 1 

c W 0 (k* ,eo,Z) (l-U sirup cosGf ’ 

sin 6 b 2 {kl,(o)V l {k* x ,ri,(o) 1 
c W o (kl,(0,Z ) (l-f7cos0) 2 t a ‘ 


(55) 


Here W a - D ] W a is dimensionless, and derivatives present in the Wronskian are now evaluated 
with respect to the normalized distance 77 . Factor ( Dj / 1 0 ) present in the 2D solution implies 
that the pressure field is evaluated per unit flow thickness t 0 in x 2 direction. Since the 
Wronskian in equation (23a) is independent of normal distance to the wall, computationally it is 
advantageous if it were evaluated at infinity where V 2 may be eliminated from the solution 
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( 56 ) 


b 2 (k*,co) 

w 0 (k;,co,z) 


J D jX oJi 


iDjXJV^k* ,ri,(o) + V'ik'l ,ri,co ) 


1 7] — >°< 


Sample Computations 

We start with a rigid boundary condition on the wall (Z — » ») and a set of values for the 
observer location and flow condition as (0 = n 12,0 - n /4,St 0 = 0.30, Uj lcj= 0.90, T R - 3.0) . 
Two linearly independent solutions V x and V 2 to equation (23a) are shown in Figures 7 and 8 . 
The first solution V x is real, and possesses a relatively larger amplitude. The second solution V 2 
is complex, and is solved subject to specified boundary conditions at 7] = 0 and °° . The missing 
initial boundary condition is solved as V 2 (0) = (1. 15751, 0.06114/). At large distance r\ from 
the boundary, both Vj and V 2 exhibit an oscillatory behavior with a regular period equal to 
2 n/(-DjXJ. 

Constant b 2 is determined from Eq. (29) as 77 — » <». Figure 9 shows that 
b 2 =(-0.18368,-0.13127/). The real and imaginary components of G N (j],(o) are presented in 
Fig. 10. Since the far- field pressure is evaluated according to equation (44) as the convolution of 
the source correlation function and the GF, the importance of the GF appears only where the 
source density is not insignificant. This, most likely, is limited to the initial several diameters 
normal to the wall where jet/wall turbulence and/or heat related sources are present. Figure 1 1 
displays an expanded view of the GF within the first 2 diameters normal to the wall. The 
amplitude is relatively larger at 77 = 1 compared to regions within the TBL, and the phase is 

constant at 0.623 Rad. 


V1O7) 



Figure 7. Solution 1/(7]) (solid line) and its derivative (dashed line) subject to the mean velocity and 
temperature profiles of Figs 5 and 6 at (0 = n / 2 , 6 = n / 4, St o = 0.30, U j / = 0.90, T R = 3.0). 
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V 2 (v) 



Figure 8. Real (solid line) and imaginary (dashed line) components of solution V 2 (r/) subject to the 
mean flow profiles of Figs 5 and 6 at (0 = n/2, 6 = 7T/4, SY o =0.30, U j Ic^— 0.90, T R = 3.0). 


V 2 (7]) e iT]D ;X~ 



Figure 9. Real (solid line) and imaginary (dashed line) components of constant 
parameter b 2 as 7] — > °o . 


Gn 



Figure 10. Real (solid line) and imaginary (dashed line) components of the Green’s function 

G N (jj,(o) at ty = n/2,e = n/4,St o =030,Uj/c„=0.90,T R = 3.0). 
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Figure 11. Magnitude (solid line) and phase (dashed line) of the Green’s function G N (r\,co) 
at (0 = nil, 6 = n/4, St= 0.30, U j lc„= 0.90,T R = 3.0). 

Next, we increase the frequency to St a = 1 .0 while the remaining parameters are kept the same as 
before. Noticeable reduction in the magnitude of the GF (see Fig. 12) indicates that sound 
generation within the TBL should be minimal relative to the previous example at Sf o =0.30 
provided that the source densities remain comparable. 

Gn 



Figure 12. Magnitude (solid line) and phase (dashed line) of the GF at 

St 0 = 1.0- with the remaining parameters the same as in Fig. 11. 


In the third example, temperature effect is examined by considering an unheated jet at T R = 1 .0 , 
while the remaining parameters are kept the same as in Fig. 1 1 . The temperature profile for this 
flow is seen in Fig. 13, and the corresponding GF solution is presented in Fig. 14. The 
magnitude of the GF at and near the wall (77 < 1.1 ) is increased considerably (nearly four times) 
relative to the heated jet shown in Fig. 11. It is noted that according to equation (55), that 
temperature effect enters the computations through factors (c m / c) as well as solution V t (rj). 
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Figure 13. Mean flow profiles for temperature 
(solid line) and sound speed (dashed line) in the 
unheated jet at T R — 1.0. 


Gn 



Figure 14. Magnitude (solid line) and phase 
(dashed line) of the Green’s function G N {r\,(ti) at 

(0 = nil, e = n/4, St o = 0.30, U j !cj= 0.90 ,T R = 1 .0). 


The effect of the observer location on the GF is examined by considering a polar angle of 
0 -n / 2 while the remaining parameters are kept the same as in example 1. Figure 15 shows an 
increase in the amplitude of the GF in a direction normal to the plate. In particular, there is 
substantial increase in the magnitude of the GF within the BL at rj < 0.13 . 


IGJ 



Figure 15. Magnitude of the Green’s function G N at 
6=71/ 2 (solid line), and 6=71 / 4 (dashed 
line), at conditions of Fig. 11. 
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The effect of wall conditions may be examined by selecting numerical values for the normalized 
impedance function Z . The real and imaginary components of the surface impedance depend 
not only on the wall conditions, but also on the external flow characteristics such as Mach 
number and sound frequency and amplitude [Refs. 24, 25]. Aside from the solid boundary 
conditions (Z — » °°) examined earlier, here we choose a pair of values for the specific resistance 
and reactance as Z = (l±0.50/). Computational results are shown in figure 16 at unheated 
conditions T R = 1 .0 , and at a polar angle of 6 =n / 4 . The amplitude of the GF is affected slightly, 
and only within the first diameter normal to the wall. The phase (shown in radians) remains 
constant for the rigid wall, but becomes a function of the source location when the impedance 
function consists of a reactive component. 


I G h \ 



1.1 - 



0.9 - 


0.8 - 

0 0.5 1 


-i i — i h 

1.5 2 


Figure 16. The effect of wall impedance on the magnitude 
(top) and phase (bottom) of the GF - (rigid: solid line); 

(1 +0.50 i : dashed line); and (1 - 0.50 i : dash-dot) at 


(0 = nil, 9 = n/4, St o = 0.30, U . /c_= 0.90, T R = 1 .0) . 


7.0 Summary 

Prediction of the aircraft noise due to the interaction of the jet exhaust with a nearby solid 
surface is of increasing interest in the design of future civil transport. Concepts such as hybrid 
wing body aircraft, over the wing engine mount, or distributed propulsion are likely to contribute 
to jet surface interaction noise. While nearby surfaces could provide significant noise shielding 
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in certain directions, they also increase noise in other directions. In this study, we presented a 
formal solution to the propagation Green’s function applicable to a rectangular jet exhaust in the 
proximity of a flat surface. The governing equation is the compressible Pridmore-Brown 
equation. Both frictional heat generation and thermal gradients in the boundary layer are 
potential contributors to noise generation and propagation. The general expression for the far- 
field noise (eq. 42) resembles that in the absence of a nearby surface - the major difference is 
concealed in the propagator, which should now be evaluated subject to prescribed boundary 
conditions on the surface. It was pointed out that aside from the scrubbing noise, scattered noise 
(or trailing edge noise) is also present on both side of the surface. A comprehensive prediction 
model needs to account for both noise components. 

Numerical results were presented that show a parametric study of the propagation Green’s 
function at selective values of the source frequency, jet temperature, observer angle, and wall 
impedance function using analytical representations of the mean flow profiles. For round jets, a 
conformal mapping to a rectangular strip needs to be carried out. In practice, a Reynolds- 
Averaged Navier-Stokes solution to the nozzle flow in the proximity of the surface (with an 
appropriate turbulence model) provides the mean flow and turbulence information required for 
modeling both momentum-flux and enthalpy-flux source components. The propagation equation 
should be solved numerically at each pair of observer angles ( 0 , </>), and at each frequency co 
subject to the local mean velocity and density profiles. In a 2D formulation, the GF is evaluated 
according to equation (51) while the mean flow is considered as a superposition of parallel slices 
in jCj direction and infinitely long in x 2 direction. When the locally parallel flow approximation 
is extended to x 2 direction the GF is evaluated according to equation (39). Here, the jet may be 
sub-divided into elements in the span-wise x 2 direction in order to account for a finite wetted 
domain with an eventual decay of the source density in this direction. 


Nomenclature 

b 2 Far field amplitude (eq. 29) 

c Sound speed 

Dj Jet Diameter 

8 0 Boundary layer thickness 

77 Normalized distance y 3 / Dj 

G Green’s function 

G, G Transformed Green’s function 

G n Normalized Green’s function 

h ' Enthalpy fluctuations 

h Enthalpy 

h' o moving frame stagnation enthalpy ( h' + V- v'j I 2) 
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V 

7 

r 

k 

K 

K o 

M 


co 

n 

P' 

P(k t ,co) 

q(y,H,T) 

p 

R 


l 

¥ 

t 

T 

1 R 

T 

u 


Vi 

e 

0 

x 

y 

z 


Normalized distance from the wall 

Specific heats ratio 

Source 

Wave number 

Wave number vector (k x ,k 2 ) 

Wave number magnitude ( CO lc ^ ) 

Mach number (U / c) 

Radian frequency (2 71 f) 

Normalized pressure fluctuation 
Acoustic pressure 

Wall-pressure wave number- frequency spectrum 

Source correlation function 
Density 

Distance \x— y\ 

Spatial separation vector 
Wall function (eq. 21b) 

Time 

Stagnation temperature ratio 

Time delay 

Mean axial velocity 

Peak velocity in the boundary layer 

Turbulent velocity component 
Mass-averaged velocity component 
Polar angle 
Azimuthal angle 

Wronskian 

Rectangular coordinates 
Source location 

Surface impedance 


Subscripts 

0 At the surface variable 

00 Ambient condition 


Superscripts 

Time average 

A Fourier transformed variable 

~ Favre average (q = pq / p) 
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Appendix-A Influence of the Solid Surfaces 

The Navier-Stokes equations when linearized about a non-radiating base flow may be written as 
a set of five linear equations [Ref. 15] 

L Vfl u fl = f v (x,t), (/t,v = 1,2,..., 5) (Al) 

where u denotes the five dependent variables 


u ll =(m l ,p',p') 
m i = pv', i = 1,2,3 


(A2) 


(Latin indices repeat from 1 to 3, and Greek indices repeat from 1 to 5). The first three equations 
in (Al) correspond to the linearized momentum equation, followed by the energy and mass 
equations in the fourth and fifth place. The Green’s Function (GF) to set of equations (Al) with 
a non-zero delta function source placed in the ath equation is denoted as g 

L vfl g fia (x,t\ y, t) = S v<J d(x-y)S(t - r). (A3) 


The source term appearing on the right hand side of (Al) is 


f v ^, t ) = ^< J +S v4 (r-l) e "^-, (i,j = 1,2, 3), v = 1,...,5 


(A4) 


<=-(pvV-pvV), 


v' 4 =(y-l)(h' + %' 2 ), v'=0. 


(A5) 


The first three source components relate to difference between momentum flux terms and their 
Favre-averaged values, while the fourth term is related to the enthalpy flux and its Favre- 
averaged value, and the fifth source term (mass equation) is obviously zero. 


Note that both momentum variable m ( . and pressure variable p' are non-linear, however, in 
practice we set m i - pv ' , and 



(A6) 


implies that p' - // in the acoustic domain where turbulence is zero. 

The field variables of interest are obtained from equations (Al) and (A3) using the source/GF 
convolution integral 
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(A7) 


Upix, t) = J J g^,(x,t I y,z)f v (y,z)dydz, 

T y 


For example, when p = A the far-field acoustic pressure is evaluated as 


p\x,t ) « p'(x,t) = J J g 4v (x,t I y,z)f v (y,z)dydz . 

T y 


(A8) 


Considering that/ 5 =0, the latter operation requires solving (A3) four times, each time by 
placing a delta function in the one of the first four equations. This process would also solve for 
other GF components that are not necessary for the purpose of evaluating p ' . It is 
computationally advantageous if the field variables were evaluated from the set of the adjoint 
equations that are governed by 

4 ? \y,z) = 8 va 8(x-y)8(t-z). (A9) 


The adjoint operator (denoted by superscript a) is readily evaluated once we multiply equation 
(Al) by g ( v a J (x,t I y, z ) , and rearrange its left hand side to appear as 


gvJ (x,t I y,z)L Vfl u^(x,t) = u M (x,t)L ( “l g ( v a J(x,t I y,z) + ) + ^F ja . 


dxj 


(A10) 


The last two terms on the RHS of (A 10) represent the bi-linear form of u and g ( °J , and 


H<7 


F ja (x,t\y,z) = v j (xd)u^xJ)g^(x,t\y,z) + u 4 g ( “J 



(All) 


We solve for u G (y,z) by placing (A9) into (A10), and integrating this equation within a large 
four-dimensional space-time volume (x,t) surrounding the source 


J J g ( v a J(x,t I y,z)f v (x,t)dxdt 


u a (y,z)+j\ (^ 4 )+ F ja . 

t x V J J 


\dxdt . 


(A 12) 


Following the divergence theorem, the volume integral on the RHS of (A 12) may be converted 
into a surface integral. When there are no solid surfaces to interfere with the sound, this integral 
is evaluated on the surface of a large 4D sphere in (x,t) domain - resulting in a null contribution 
due to the vanishing of the field variables u on such a surface. Consequently (A12) becomes 


u M (x,t) = J J gv/liy,* I x,t)f v (y,z)dydz . (No solid surfaces) (A 13) 

T y 


Equations (A7) and (A13) point to the reciprocity of the GF in the absence of solid surfaces, i.e., 
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gyjiyj I X,t) = g^ v (x,t I y, r) . 


(A 14) 


By setting ji = 4 we evaluate the far- field pressure. It is seen that the advantage of (A13) over 
(A7) is that we now place a source only in the fourth equation within set (A9) to compute the 
required GF components for pressure. Equation (A 13) represents the fundamental solution to the 
governing equations. In the presence of solid boundaries in a finite domain, the general solution, 
when formulated in terms of the adjoint GF, is complemented with surface integrals as shown in 
(A12). The reciprocity condition (A14) would not necessarily be satisfied when reflecting 
boundaries are present in the acoustic medium unless surface integrals are forced to disappear by 
imposing conditions on the adjoint GF. 


To investigate the surface integrals further, we allow volume x(t), externally and/or internally, 
be bounded by surface s(t) with a unit outward normal n (Fig. Al). The first term within the 
bracket on the RHS of (A12) is expressed as the sum of two terms following the Leibniz’s rule 



y J u M gvo (x,t\y,t)dx- j n- V s (. x,t\y , r)ds 

dt m s(t) 


(A 15) 



Figure Al. Solid surfaces within the acoustic medium 


where V s (x,t) is the velocity at an arbitrary point on the surface s(t ) . We now integrate both 
sides of (A15) within the time interval [-T, T] for some large time T. The first terms on the RHS 
of (A 15) vanishes due to the initial condition (in a remote past time), and the causality condition. 
Inserting (A15) into (A12) and applying the divergence theorem to the last term in (A12), we 
find (upon interchanging dummy variables) 

T 

u M (x,t) = jjg < ^(y,T\x,t)f v (y,r)dydr + j dr J n j \y.u v g^{s,r\x,t)-F jll {s,r\x,t)\ds (A16) 

T y -T s(t) 

Integral equation (A 16) represents a formal solution to the problem. The two integrals represent 
contributions to the acoustic field due to the volume and surface sources, respectively. In 
particular, the acoustic pressure is solved when /i = 4 . Moving the partial derivatives from 
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source to the GF and/or imposing boundary conditions on the surface integrals may achieve 
further simplifications. For example, placing f v from (A4) into (A16), and using the identity 

gjf(y,r I x,t)^-e"j = -^-(g ( “J(y,T I x,t)e")~ e' v 'j ^ g ( v a J (y ,T I x,t) 
followed by application of the divergence theorem on the first term on the right shows that 


“„<*.<) =-jj 


< 3g M 


T y 
T 


dv„ 


-^-(r-D^gS 

dyj dyj , 


1 

e"(y,r)dydT+ J drjn j g ( ^(s,T\x,t)e"jdsdT 


-T s 


(A 17) 


+ J drjnj[VjU v g ( v a M \s,T\x,t)-F^(s,T\x,t)\ds. 


-T s 


Further, substituting for Fj from (All) into this last expression and imposing the no-slip 
boundary condition on a hard surface 


(v;-v.).«. = o, 


(A 18) 


results in the following integral equation for pressure 


1 

Pe(x,t ) = - J J nyobr I x,t)e" J (y,t)dydz+ J dr J tijg^is, x I x,t)e"jds 

T y -T s 

T 

- J dT\n j [_p' e {s,T)g^{s,T\x,t) + m j {s,T){g^ + c 1 g^)\ds 


(A 19) 


-T s 


where 




p(«) — 




dy. 


(A20) 
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Appendix-B Vanishing of the Surface Integrals 

Here we show that the boundary condition (21) eliminates the surface integrals in the Green’s 
function formulation of the acoustic field for bounded media. We start with equation (1), and 
apply Fourier transform with respect to variables xi, X 2 and t as defined in eq. (7) with the over 
hat carrot denoting a transformed variable 

d 2 tf ~ - 

j^+f(k t ,x 3 ,G))# = A, (Bl) 

where function / was defined in eq. (1 1) and 

d(k,x 3 ,co) = n'(k ,x 3 ,co) , A (k,x 3 ,co) = -r . (B2) 

-(0 + k x U{x 3 ) c(x 3 )(-a> + k x U (x 3 )) 


Multiply eq. (9) by & , and equation (Bl) by G , subtract the two expressions, and integrate the 
result with respect to x 3 



d 2 tf 

dx 2 


-i? 


d 2 G 

dx 2 j 


dx 3 = J^GA-t?5(x 3 -y 3 )jdx 3 , 

o 


(B3) 


or 

= d(k,y 3 ,co)- ^G(k t ,x 3 ,y 3 ,co)A(k,x 3 ,Q))dx 3 . (B4) 

o 

Upon recognizing that dd / dx 3 = and rJG / dx 3 = —iXJ& at infinity, the upper limit of the 

expression on the left hand side of (B4) vanishes. On the right, since the GF is self-adjoint we 
switch x 3 and y 3 in the argument of G , and subsequently we substitute expressions for A from 

(B2), and G from (10) into (B4) to show that in view of (19), the right hand side of this equation 
also vanishes. 




dx, dx 


3 J 


Consequently we arrive at 


g|^— ^|^ = 0, x 3 = 0 . 

dx 3 dx 3 

Using surface boundary condition (18) in (B5) shows that 


y/G - 0 , x 3 -0 

dx 3 

where function y/{k x ,(0,Z) is given in (21b). 


(B5) 


(B6) 
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